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Resumen 



La dinamica de un gas autogravitante de neutrones en presencia de un 
campo magnetico, es estudiada tomando como ecuacion de estado del gas de 
neutrones la obtenida en un estudio previo [2]. Se trabaja con un espacio 
Bianchi I caracterizado por una metrica tipo Kasner que permite tomar en 
cuenta la anisotropia que introduce el campo magnetico. El sistema de ecua- 
ciones de campo de Einstein Maxwell para este gas se convierte en un sistema 
de ecuaciones dinamicas en un espacio de fase en 4D. Se obtienen soluciones 
numericas de dicho sistema. En particular se encuentra una solucion singu- 
lar tipo punto para diferentes condiciones iniciales. Fisicamente esta solucion 
singular podra asociarse al colapso de un volumen local de materia neutronica 
en el interior de una Estrella de Neutrones. 



Abstract 

The dynamics of a self-gravitating neutron gas in presence of a magnetic 
field is being studied taking the equation of state of a magnetized neutron 
gas obtained in a previous study [2]. We work in a Bianchi I spacetime 
characterized by a Kasner metric, this metric allow us to take into account the 
anisotropy that introduces the magnetic field. The set of Einstein-Maxwell 
field equations for this gas becomes a dynamical system in a 4-dimensional 
phase space. We get numerical solutions of the system. In particular there 
is a unique point like solution for different initial conditions. Physically this 
singular solution may be associated with the collapse of a local volume of 
neutron material within a neutron star. 
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0.1 Introduccion 



Entre los objetos que habitan en el universo uno de los mas interesantes y 
que poseen mayor riqueza de fenomenos fisicos son las Estrellas de Neutrones 
(ENs). El estudio de este tipo de objetos presenta una gran variedad de 
atractivos para los fisicos, tanto teoricos como experimentales. 

En las ENs concurren una amplia gama de condiciones fisicas extremas 
como son: altas densidades, elevados campos magneticos y el hecho de ser 
una de las cinco configuraciones estelares (Enanas Blancas, Estrellas de Neu- 
trones, Huecos Negros, Estrellas Supermasivas y grupos de estrellas rela- 
tivistas) en las cuales es preciso tener en cuenta los efectos de la Relatividad 
General (RG) p. Su complejidad hace que ademas de la Gravitacion haya 
que tener en cuenta las otras tres interacciones fundamentales (fuerte, debil 
y electromagnetica). 

En el alio 1934 se predijo por primera vez, basandose en calculos teoricos, 
la existencia de una ENs y a partir de esa fecha se estudiaron muchos mod- 
elos para describirlas, pero no es hasta el afio 1967 en que se verifica obser- 
vacionalmente por azar la existencia de un objeto con las caracteristicas de 
una ENs. 

Las ENs no se habian observado antes debido a que no emiten en el 
visible, por lo cual fue necesario que pasaran mas de 30 afios y que la tec- 
nologfa permitiera la construccion de radiotelescopios capaces de detectar las 
primeras ENs. El descubrimiento de la primera ENs provoco un aumento en 
el mimero de trabajos teoricos encaminados a describir los procesos fisicos 
que ocurren en el interior de estos objetos. 

Desde hace veinte afios los estudios de las Estrellas de Neutrones se nu- 
tren de cuantiosas observaciones que se obtienen gracias a observatorios que 
detectan emisiones de rayos X y radiacion gamma colocados en satelites. 
Entre los observatorios mas conocidos estan el Telescopio Espacial Hubble, 
el Observatorio de Rayos X Chandra y el Observatorio de Rayos Gamma 
Compton. 

Actualmente las principales interrogantes concernientes a la fisica de las 
ENs son: 

• Determinacion de cotas para la masa y el radio. Determinados por 
ecuaciones de estado mas realistas. 

• Estructura interna de las Estrellas de Neutrones y sus consecuencias 
observacionales: estudio de los Sideramotos o disminucion del perfodo 
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de rotacion que pudiera deberse al cheque entre la corteza y liquido 
superfluido. 

• Origen, evolucion y cotas de los campos magneticos que ellas soportan. 

• Determinacion de las propiedades de la materia en regimenes criticos, 
alta densidad por encima de la densidad nuclear. 

• Existencia de ondas gravitacionales como consecuencia de los marcados 
efectos relativistas que traen consigo las ENs. 

Al tratar de dar respuesta a estas preguntas se tropiezan con muchas di- 
ficultades por lo que es preciso hacer modelos que describan de forma aprox- 
imada los procesos que ocurren en el interior de las ENs. Muchos de estos 
modelos se basan en tomar en cuenta diferentes ecuaciones de estado y luego 
comparar los resultados con los que se obtienen por medio de las observa- 
ciones. 

Como los campos magneticos en las ENs son elevados resulta razonable 
estudiar las ecuaciones de estado del gas de neutrones magnetizado. En 
[2] se ha estudiado el gas de neutrones relativista en presencia de campo 
magnetico interactuando con el a traves del momento magnetico anomalo. 
Uno de los resultados obtenidos en este estudio es que las presiones presentan 
anisotropfas y el gas colapsa en la direccion del campo. 

Debido a que los sistemas con anisotropfas, como el descrito en [2] en 
general no son estables en el universo un objeto de este tipo deberfa colapsar 
o realizar una transicion a un estado estable. 

Incluir la gravedad nos permite obtener informacion sobre el papel que 
ella juega cuando aparecen estos tipo de anisotropfas en las presiones. Esta 
serfa la principal motivacion de este trabajo y nos llevarfa a responder las 
siguientes preguntas: 

• iLa. Gravedad frenara o estimulara el colapso? 

• iPodia este sistema evolucionar a una configuracion estable? 

Con el objetivo de responder estas interrogantes, estudiaremos en esta tesis, 
el gas de neutrones magnetizado autogravitante. Utilizaremos una metrica 
no estacionaria para obtener la evolucion en el tiempo del sistema. 

Por supuesto la formulacion de un modelo de una estrella neutronica en 
presencia de campo magnetico requiere del uso de una metrica mucho mas 
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complicada, que se traduce en un problema niimerico mas engorroso. Nuestro 
modelo aunque simplificado puede dar informacion cualitativa muy intere- 
sante sobre los procesos dinamicos que le ocurren a un volumen local en el 
interior de una ENs y constituye un primer paso en el estudio de un sistema 
tan complejo. 

En un trabajo anterior [3] fue estudiado el gas de electrones magnetizado 
utilizando igualmente un espacio-tiempo Bianchi I y una metrica Kasner. 
Este trabajo pretende entonces hacer un estudio similar para el gas de neu- 
trones. 

La tesis tendra 2 capftulos introductorios que consideramos basicos para 
entender la problematica que nos ocupa. los capftulos 3 y 4 se dedicaran a 
la parte original de la tesis. 

El contenido esta distribuido de la siguiente forma: 

Capi'tulo 1 se describen las caracten'sticas generales de las ENs asi como 
su formacion y evolucion, ademas se discute el rol de los campos magneticos 
en estos sistemas y se discuten algunas evidencias observacionales que existen 
actualmente sobre las ENs. 

Capitulo 2 se obtienen las ecuaciones de estado de un gas degenerado de 
neutrones en presencia de un campo magnetico interactuando con el a traves 
del momento magnetico anomalo de los neutrones. 

Capitulo 3 se utiliza la RG, introduciendo variables dinamicas para 
obtener un sistema de ecuaciones diferenciales que describen la evolucion 
de un volumen local de una ENs. 

Capitulo 4 se estudian las soluciones del sistema de ecuaciones en funcion 
de diferentes condiciones iniciales asi como las trayectorias en el espacio de 
fase. Posteriormente se exponen las conclusiones. 
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Capitulo 1 

Estrellas de Neutrones. 



1.1 Caracterfsticas generales. 

El fisico sovietico L. D. Landau fue el primero en sugerir que cuando la 
densidad de la materia es muy grande, los electrones se ven obligados a 
fusionarse con los protones. Landau predijo que se llegan'a a una nueva 
configuracion de equilibrio, en la que la densidad de la materia es tan alta 
que los micleos atomicos en contacto forman un gigantesco niicleo. 

En el afio 1932 el neutron fue descubierto por Sir James Chadwick lo cual 
aclaro el problema de la evolucion estelar a altas densidades. La existencia 
de un objeto estelar compuesto fundamentalmente por neutrones fue pro- 
puesto de forma teorica por primera vez en 1934 por Walter Baade y Fritz 
Zwicky |1] en un articulo sobre la naturaleza de las supernovas, sugiriendo 
que las ENs se formaban en este tipo de explosiones. Son de seiialar los tra- 
bajos de Oppenheimer y Volkoff realizados en 1939 [S], en los cuales estudian 
las posibles configuraciones de equilibrio de una ENs utilizando un modelo 
de estrella en el cual la presion del gas degenerado de neutrones compensa 
la presion gravitacional como mismo ocurre con el gas de electrones en las 
Enanas Blancas (EEs). 

Actualmente se consideran como Estrellas de Neutrones (ENs) [B] a las 
estrellas con masas del orden de 1.5 veces la masa del sol (M = IjSM©), 
radios de alrededor de 12 km [R ~ 12 km) y densidades comparables con 
las del nucleo atomico (p ~ lO^^g/cm^). La temperatura superficial de una 
estrella de neutrones se encuentra entre 3 x 10^ K y 10^ K y la temperatura 
central es del orden de 10^ K. La energia de Fermi (fip) la podemos estimar 
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conociendo que: 



iVm„ N p 6.022 X lO^^m-^ para p,„„ ~ 10^° kg/rn^, 

' V V rrin [4,21 X para p^ent ~ 7 x 10^^ kg/m^, 

donde p es la densidad, es el niimero de particulas, m„ es la masa del 
neutron y es el volumen. La energi'a de Fermi es [7j: 

, oN./n/nx I 1.06 X 10"^^ J, en la superficie, 
2mn V [1,80 X 10-1^ J, enelcentro, 

y la temperatura de Fermi (Tp = ^) queda: 

kB 

7.69 X 10^ K, en la superficie, 
1,30 X 10^2 K, enelcentro. 

Para que un gas fermionico se pueda considerar degenerado tiene que cumplirse 
T 

la condicion — ^ 1, como vemos en la ENs esto se cumple para regiones 
Tp 

T 

cercanas al centro donde — ~ 10~^. 

Aunque la mayor parte de la estrella esta formada por neutrones, existen 
ademas protones y suficientes electrones para garantizar la neutralidad de la 
estrella. Producto de las altas densidades tambien se pueden encontrar toda 
una amplia gama de particulas en su interior. La estructura interna de una 
ENs consta de una serie de capas de diferentes composiciones, la misma se 
divide en 5 partes fundamentalmente [B]. Una represent acion esquematica se 
muestra en la figura fll.ip . 

La atmosfera y la envoltura contienen una cantidad despreciable de masa, 
pero la atmosfera juega un papel importante en la forma del espectro de los 
fotones que salen de la estrella, y la envoltura tiene una influencia crucial 
en la liberacion y el transporte de la energi'a termica de la superficie de la 
estrella. La corteza se extiende entre 1 a 2 km por debajo de la superficie 
y esta formada fundamentalmente por nucleos atomicos. El nucleo de la 
estrella contiene el 99% de la masa de la estrella. El nucleo exterior consiste 
en una combinacion de nucleones, electrones y muones. En el nucleo interior 
pueden existir particulas exoticas, se piensa que incluso puede ser posible la 
transicion a una fase mixta de materia hadronica y quarks libres. 
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1.2 Formacion y evolucion. 

Cuando las estrellas llegan al final de su vida, producto del agotamiento 
del combustible que mantiene las reacciones termonucleares necesarias para 
garantizar el equilibrio, evolucionan hacia un estado final dependiendo de su 
masa. Algunas explotan dispersandose en el medio estelar, otras se contraen 
formando enanas blancas, ENs o agujeros negros. 
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Las estrellas de neutrones surgen producto del colapso gravitatorio del 
niicleo de una estrella masiva (> SM©). Al final de su vida una estrella 
masiva expulsa la mayor parte de su masa en una explosion conocida como 
supernova tipo II, es precisamente en esta explosion donde surgen las ENs. 
En una estrella ocurren reacciones de fusion en las cuales el hidrogeno se 
transforma en helio, este por ser mas pesado se va acumulando en el nucleo, 
cuando este alcanza una temperatura de 10^ K los 3 procesos alfa comienzan, 
formandose g^C mediante las reacciones: 

^He+^He ^ ^Be +^ He ^ + ^ + 7.367MeV. 

Reacciones como estas se repiten formandose elementos mas pesados (^'^O, 
^^Si). Mientras mas pesados sean los elementos producidos en cada reaccion 
menor es la energia que se libera, las reacciones continiian hasta que se forma 
^^Fe que es el nucleo con menor energia por nucleon. La energia no puede 
continuar produciendose por la fusion y el nucleo de la estrella formado por 
hierro crece. Este aumento contimia hasta sobrepasar el valor de 1,4Mq 
(limite de Chandrasekhar) , en este punto el niicleo comienza a colapsar, pues 
la presion del gas degenerado de electrones no es suficiente para contrarrestar 
la fuerza gravitatoria y las capas superficiales son expulsadas hacia afuera, 
esta explosion es la supernova. 

En el nucleo de la estrella, protones y electrones son comprimidos formandose 
neutrones (n) y neutrinos (z/g)- Producto del colapso la densidad aumenta 
y los micleos atomicos se vuelven mas ricos en neutrones, los cuales comien- 
zan a escapar de los micleos. Para densidades superiores (p ~ lO^^g/cm^) 
los micleos remanentes se desintegran quedando la materia compuesta,en 
su mayor parte, por neutrones. Los neutrones libres son inestables con un 
tiempo de vida media cercano a los 15 minutos decayendo a traves de la 
reaccion /3~ : n — > p + e~ + z/~. 

En la materia a altas densidades el decaimiento I3~ no ocurre gracias al 
Principio de Exclusion de Pauli [8]: Si todos los niveles energeticos bajos de 
los electrones y protones estan ocupados, solo los neutrones suficientemente 
energeticos podran decaer. 

Por otro lado, los protones y electrones energeticos podran combinarse 
mediante la reaccion (3^ : e~ + p — > n + v^- Entre los procesos (3^ y (3~ 
habra un equilibrio si los potenciales quimicos de los neutrones, electrones 
y protones cumplen la relacion: /i„ = /ie + /ip. El destino final del niicleo 
de la estrella esta determinado por su masa. Si esta se mantiene por debajo 
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de 3 Mq se habra formado una ENs, si es superior el niicleo continuara su 
colapso formandose asi un agujero negro. 



1.3 Campos magneticos intensos. 

Las ENs son uno de los objetos astrofisicos a los cuales se les atribuyen los 
mayores valores de campos magneticos que existen en la naturaleza. Estos 
superan por muchos ordenes de magnitud a los campos que son familiares 
para los seres humanos. Por ejemplo, el campo terrestre es 0, 6 G en los polos, 
el campo permanente mas intenso que se ha obtenido en el laboratorio es de 
10^ G mientras que el campo de una ENs puede llegar hasta 10^^ G. 

Existen varias hipotesis para explicar la formacion de los campos magneticos 
en las ENs. El modo mas simple de explicar su surgimiento es producto de 
la compresion del fiujo magnetico durante el colapso del nucleo de la estrella 
progenitora: 

Ro \2 



Bens=(^) So ~ (10"- 10^3) G, 



donde Bens es el campo de la ENs, Bq es el campo del nucleo de la estrella 
progenitora, Rq = (10^ — 10^) km son los valores entre los que se encuentran 
los radios del nucleo de la estrella progenitora y Rens = 10 km es el radio 
tipico de una ENs. Con este mecanismo se explican los campos magneticos de 
los radiopulsares, sin embargo hay observaciones que arrojan valores mayores 
en los campos (~ 10^^ G) asociados a objetos Uamados Magnetars. Para 
explicar el surgimiento de tan elevados campos se han propuesto mecanismos 
mas complejos, como son los procesos tipo dinamo que consisten basicamente 
en la generacion de corrientes internas producto de la rotacion. La explicacion 
y determinacion del origen, valor y evolucion de los campos magneticos es 
un problema abierto para los astrofisicos y fisicos de parti'culas. 

La cota maxima de los campos magneticos la podemos estimar exigiendo 
que la energia magnetica sea menor que la energfa gravitatoria: 

Stt 3 R 

Estos campos magneticos se encuentran en ENs jovenes pues a medida que 
la estrella envejece pierde energia rotacional y los campos magneticos dis- 
minuyen hasta valores de 10^ G. 
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1.4 Evidencia observacional. 

En 1967 se confirma observacionalmente la existencia de las ENs k 
los esfuerzos de Joselyn Bell, estudiante de doctorado del profesor Antony 
Hewish de Cambridge, que descubre la primera fuente de radio (Pulsar), 
en el Radio Observatorio de dicha Universidad. Este primer pulsar Ueva 
el nombre de CP 9119. Actualmente se conocen mas de 1600 pulsares, 
que no son mas que ENs que emiten radiacion electromagnetica de forma 
periodica. El periodo mmimo de una ENs lo podemos estimar igualando la 
fuerza centrifuga en el ecuador con la fuerza gravitacional: 



Este valor tan bajo del periodo es uno de los principales indicadores de que 
estamos en presencia de una ENs. Los perfodos de los pulsares observados 
van desde 8,51s (PSR J2144-3933 [9]) hasta 716 Hz ^ 1.3 ms (PSR J1748- 
2446ad) [10] que es el pulsar mas rapido que se conoce. Recientemente ha 
sido descubierto una ENs con una frecuencia de 1122 ± 0.3 Hz, la misma 
fue observada en el transiente de rayos X XTE J1739-285 [11]. Objetos 
compactos que roten con perfodos menores que 0, 5 ms han hecho pensar en 
la posibilidad de la existencia de estrellas de quarks (EQs). 

Las masas observadas para las ENs se encuentran entre 1 Mq—2, 1 Mq [T2] 
con un valor promedio de 1.35 M0, asi tenemos al pulsar PSR J0751+1807 con 
una masa de 2, 1 ± 0, 2 M0, el caso mas tipico se puede encontrar en el pulsar 
binario PSR 1913+16 |I2] en el cual las masas son (1.3867 ± 0.0002) M© y 
(1.4414 ±0.0002) Mq respectivamente. La masa mas pequefia es la del pulsar 
J1756-2251 [I3| con un valor de (1.18 ± 0.02) Mq. 

Las mediciones de los radios de las ENs son muy imprecisas (9 — 15 km) 
pues no se conoce exactamente la composicion qufmica de la atmosfera ni la 
distancia exacta de las estrellas. 
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Capitulo 2 

Gas de neutrones. 



2.1 Ecuacion de Dirac para particulas neu- 
tras con momento magnetico anomalo. 

El gas de neutrones en presencia de un campo magnetico ha sido estudi- 
ado en [2], [H]. En esta seccion nos limitaremos a exponer los resultados 
fundament ales que nos seran utiles en la tesis. 

Trabajaremos en el Ensemble Gran Canonico, considerando un subsis- 
tema (en el contexto Astrofi'sico pudiera ser un elemento de volumen dentro 
de la estrella), el cual se encuentra bajo la influencia del campo magnetico 
creado por el resto del sistema H. Producto de la presencia de este campo 
el subsistema se polariza, creandose una magnetizacion en el medio (gas de 
neutrones) que satisface la relacion H = B — AnAi . El campo H es externo 
al subsistema mientras que B es externo a cualquier particula que se escoja 
dentro del subsistema, pues esta siente en adicion a H, la contribucion AnAd 
del resto de las particulas del subsistema. 

Para encontrar las ecuaciones de estado del gas de neutrones es necesario 
calcular los autovalores de la energia de las particulas que conforman el sis- 
tema. Para esto se resuelve la ecuacion de Dirac para particulas neutras con 
momento magnetico anomalo: 

(7m^m + m + iqaf,xFf,x)'ip = 0, (2.1) 

donde a^x = |(7^7a — 7a 7^*) es el tensor de spin y Fx^ es el tensor del campo 
electromagnetico. 
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En la ecuacion fl2.ip se ha tornado el convenio h = c = 1, en este sistema 
de unidades [L] = [T] = [M~^]. A lo largo de la tesis emplearemos este 
convenio en el manejo de las ecuaciones por lo que en ninguna de estas 
apareceran factores que contengan a ^ o a c. No haremos uso del convenio 
cuando estemos estimando ordenes de magnitudes, expresando estas en el SI 
o en el CGS lo cual sera evidente por sus unidades. 

De la solucion de la ecuacion (12.11) se obtienen los siguientes autovalores 
de la energi'a [H], [T5] : 



Er,ip,B,r]) = sJpl + {^/pl + ml + r]qBy, (2.2) 

donde P3,pj_ son respectivamente las componentes del momentum en las di- 
recciones paralela y perpendicular al campo magnetico (B), rrin es la masa del 
neutron, q = — 1.91/i7v es el momento magnetico del neutron (/xtv = e/2mp es 
el magneton nuclear), 77 = ±1 son los autovalores de correspondientes a las 
dos orientaciones del momento magnetico con respecto al campo magnetico. 
No hemos tenido en cuenta las correcciones radiativas por lo cual considerare- 
mos que el momento magnetico se mantiene constante para campos elevados. 
El gran potential termodinamico toma la forma 

n=-kT\nZ, (2.3) 

donde k es la constante de Boltzman, T es la temperatura, Z = Tr{p) es la 
funcion de particion del sistema j p = e~*^^~^^''/'^^, H es el Hamiltoniano, 
el potential quimico y el operador numero de particulas. 



2.2 Ecuaciones de estado. 

La forma del tensor energia-momentum de la materia en un campo magnetico 
externo constante es [16]: 

Este tensor en el caso de limite de campo cero reproduce el tensor T*^ = 
P5\ - (P + U)5\6\ de un fluido ideal. 
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Las componentes de este tensor son: 

r\ = P3 = -f^, (2.5a) 
r\ = T\ = -Q-BM, (2.5b) 
T\ = -U = -TS -fiN-n, (2.5c) 

donde S es la entropia, = —dQ/dfi es la densidad de particulas, ^A = 
—dQ/dB es la magnetizacion, U es la densidad de energia y P3 es la presion 
en la direccion del campo. 

La expresion para el potencial termodinamico consta de dos terminos: 



Vn, 



el primero es la contribucion estadistica y el segundo la del vacio [2] , explicitamente 
tenemos: 

1 ■ r°° 

^sn = -^^ Yl / P±^P^^P3Hf^if^n,0f~{f^n,0]^ (2-6) 
^ f?=l,-l 

donde ^ = l/fc^P, /^(/^n,0 = (1 +e*^^"^^"^^) representan las contribuciones 
de las particulas y de las antiparticulas. La expresion para el termino de 
vacio esta dada por la expresion: 

1 ■ r°° 

^Vn = / PA^dpA_dp2.En, (2.7) 

que es divergente pero se puede regularizar obteniendose que para B < 10^^ G 
no es importante su contribucion [2], por lo que no se tendra en cuenta en 
nuestro estudio. 

La expresion fl2.6p se puede integrar facilmente en nuestro caso (T = 0) 
obteniendose: 



vtsn — —A 

ri=l-l 



fif^ , (1 + v/3){5vP - 3)/i/, , (1 + r](3)^{3 - y/S) ^ y/Sfi' 
12 + 24 ^ 24 " 6 



(2i 

donde hemos introducido las notaciones: 



^ .s., = arcsm , // = 

2 \ fi J nin 



fr, = - (1 + r]/3y, Sr, = - - arcsm ( — — ) , // 
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In 



+ fv 



P 



B 



mZ 



siendo Be = rrin/q — 1.56 x 10^° G el campo critico y A = ^ oJl o = 4.11 x 

A partir del Potencial Termodinamico Q podemos obtener todas las can- 
tidades termodinamicas del sistema. En particular la densidad de neutrones 
y la magnetizacion, de esta forma obtenemos que = A^orAr y A/i = ^AoTM 
donde A'o = X/rUn, J^q = Nqq y los coeficientes Fjv, Tm vienen dados por: 



N 



M 



E 

'7=1,-1 



3 2 2 ^ 



■'7=1-1 



6 



g^'7 



Teniendo en cuenta lo anterior, (12.51) y (12.81) podemos escribir las ecua- 
ciones de estado para un gas relativista degenerado de neutrones en un campo 
magnetico externo en la forma: 



U 

P 
M 



fir^N + Q = Xru{f3,fi), 
BM = X/3TMi/3,fi), 



(2.10a) 
(2.10b) 
(2.10c) 



donde 

Tp = 
Tu = 



»7=1,-1 



12 + 24 ^ 24 6 
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Capitulo 3 

Relatividad General, metrica 
Kasner. 

3.1 Relatividad General 



Debido a las altas densidades que se dan en las Estrellas de Neutrones los 
efectos de la Relatividad General son muy marcados por lo que se hace im- 
prescindible introducir las ecuaciones de Einstein para la descripcion de este 
tipo de sistema. Una prueba de la accion del campo gravitatorio en este sis- 
tema la podemos obtener facilmente estimando la velocidad de escape para 
una ENs tipica (M = M©, R = 10 km): 



donde c es la velocidad de la luz y Gn es la constante gravitacional. Si la com- 
paramos con la de la tierra (11.2 kms~^) e incluso con la del sol (617.5 kms~^) 
o la de una EB (~ 0.02 c) podemos ver que la velocidad de escape de una 
ENs es mucho mayor. 

Asumiremos que nuestro sistema queda descrito por las ecuaciones de 
Einstein las cuales vinculan el contenido de energia-materia del sistema con 
la geometria del espacio-tiempo [1]: 




(3.1) 



donde k 



% = SttGn en el SI k tiene la expresion k = — - — , Gn = 6.67 x 
.^kg"^s~^, G^u = R^iv — \RgiJLv es el tensor de Einstein el cual viene 
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-11 



m' 
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determinado por el tensor de Ricci y por el escalar de Ricci R = -R^^, 
estos ultimos dependen de segundas derivadas de la metrica. 

p pa p« I p« p/3 p/3 pa (o n\ 

las cantidades F"^ son los indices de Cliristoffel, que dependen de primeras 
derivadas de la metrica por la formula, 

r^i/ = —{9f3^l,u + gup,^i - gfiu,i3)- (3.3) 

T'^ es el tensor energi'a impulso (2.4) que en nuestro caso lo podemos escribir 



1/ 

como: 



en el cual 



= (t/ + P)n^n, + P5t + n^., P = p-^^, (3.4) 



W, = dzag[U,U,-2U,0], U = -^, U^^ = 0. (3.5) 



Para hallar la forma del tensor de Einstein debemos tomar una metrica 
apropiada, como ya dijimos es conveniente en el caso que nos ocupa elegir 
una metrica tipo Kasner pues esta es compatible con la anisotropia del gas 
de neutrones magnetizado: 

ds"" = Aitfdx^ + Bitfdy^ + Citfdz^ - dt", (3.6) 

con la metrica (13.61) y el tensor (13. 4p se obtiene que las componentes no nulas 
de la ecuacion (13. ip son: 

-G^. = |§ + | + § = -^(^^-^-^)' (3.7a) 
AC A C 

-G\ = ^ + I + ^ = -^(^^-^-^)' (3.7b) 

AB A B , , 

„, Ab Ac bc 

La notacion del punto significa derivada respecto el tiempo, por ejemplo 
• dA ■■ d^A 

A — A — 

dt' dt^' 
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La ecuacion de Einstein implica la conservacion de la energfa: 

r^''.^ = 0. (3i 

El punto y coma denota la derivada covariante, por ejemplo: 



q-t^ — _ -pi 7"^ _|_ -pt^ q-i q^ 



Tomando en cuenta la ecuacion (13. 8p obtenemos: 

i,.£(!^-(^ + |)(-BA1+p + C), (3,0) 

To memos ademas las ecuaciones de Maxwell: 

F'^^, = 0, F[^.;«]=0, (3.11) 

asi obtenemos: 

Hemos obtenido las ecuaciones (13.71 13.101 13.121) en las que tenemos derivadas 
de primero y segundo orden de las funciones B^C y U . Para obtener un sis- 
tema de ecuaciones diferenciales de primer orden cuya integracion numerica 
resulte mas sencilla introduciremos un conjunto de nuevas variables. 



3.2 Variables dinamicas. 

Para la descripcion dinamica de un elemento de volumen utilizaremos las vari- 
ables que usualmente se usan en cosmologica para caracterizar la evolucion 
del universo. Asi tomaremos un sistema comovil en el cual la cuatrivelocidad 
es = (5f y cumple que = — 1. 

Introduciremos las siguientes variables definidas a partir del vector u°' 

hap = QalS + UaUp, (3.13) 

haf3 se denomina tensor de proyeccion si este se contrae con otro tensor, 
como resultado lo proyecta sobre un 3-espacio ortogonal a la 4-velocidad m", 
cumpliendose que 

KK = K, /iA/3 = 0, h\ = ?,. (3.14) 
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Tambien se puede descomponer la derivada covariante Ua-^p en las siguientes 
partes irreducibles, 

= + + ^0/iQ/3 - (3.15) 

donde cTq/j es un tensor simetrico de traza nula, ujai3 un tensor antisimetrico, 

y M^cr^/j = u'^ujafj = 0, siendo, 

e = (3.16a) 

tta = Ua;/3U^, (3.16b) 
0-^/3 = ^{Ua;i,h''p + Up.^h^'^) -^Qha0, (3.16c) 

Cfa/S = ^{Ua;,ih^(3-Up.f,hf'J. (3.16d) 

La cantidad O se llama escalar de expansion, la misma caracteriza la veloci- 
dad a la cual el elemento de volumen ortogonal a u'^ se expande o contrae, 
tta es la 4-aceleraci6n del fiuido. El tensor aab es el rango de deformacion, 
el cual describe la manera en la que el elemento de volumen ortogonal a la 
cuatrivelocidad (u^) cambia su forma, la magnitud de este tensor viene dada 
por: 

= ^a^fsa''^, (3.17) 

y u!ai3 es el tensor de vorticidad, el cual es una medida de la rotacion presente 
en la materia y es conveniente definirlo vectorialmente como: 

= ^^^''y'u^uj.s, (3.18) 

aqui 7] es el tensor totalmente antisimetrico, ademas se satisface que u'^Ua = 
0, y la magnitud de la vorticidad viene dada por, 

^2 = ^a^^ = (3.19) 

si la vorticidad es cero entonces se dice que el vector del campo es irro- 
tacional. Esto es algunas veces conveniente para combinar la magnitud del 
rango de la deformacion con la expansion, quedando el tensor de expansion 
como, 

= (Tap + -QKf). (3.20) 
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En el caso de la metrica fl3.6l) y la 4-velocidad = 5", la 4-aceleraci6n se 
anula, el escalar de expansion 6 y el tensor de deformacion toman la 
forma: 

= ^ + | + § , a^,=diag\al,alal,% (3.21) 



donde 



2A B C 
~ 35 ~ SC' 



(3.22a) 



Haremos el siguiente cambio de notacion para mayor comodidad en el 
tratamiento de las ecuaciones: 

a: -.Si, a^^^S2, a^^Sa. (3.23) 

Haciendo uso de estas variables dinamicas y teniendo en cuenta las ecuaciones 
f l3.10p y f l3.12p obtenemos el siguiente sistema de ecuaciones: 

(t/ + p-^i3A^)e-SA^S3, (3.24a) 
^ - eS„ (3.24b) 



u = 




E2 = 






2 


S3 = 


—1 




3 


e = 


«:( 




2 


/3 = 


3' 



M - 0S3, (3.24c) 

k{BM + ^{U -p)) -Q\ (3.24d) 

S3 -29). (3.24e) 
Ademas de (3.7d) se obtiene que 

92 

- S^ - S2S3 + — -J:1 = kU. (3.25) 

Notemos que hemos obtenido un sistema de ecuaciones diferenciales no 
lineales de primer orden en las variables U, (3, 6, S2, S3 mas la ligadura fl3.25p . 
La solucion de este sistema describe la evolucion dinamica de un volumen 
local de un gas de neutrones magnetizado que podn'a ser un volumen local 
en el interior de una ENs. 
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Capitulo 4 

Ecuaciones dinamicas. 



4.1 Variables adimensionales. 

En este capitulo estudiaremos la dinamica del gas de neutrones magnetizado 
que equivale al volumen local de una ENs. Nuestro interes fundamental es 
responder la pregunta inicial de si en el marco de la Relatividad General el 
colapso obtenido en [2] se mantiene. 

Para realizar la descripcion dinamica del volumen local de la estrella de 
neutrones introduciremos un conjunto de variables adimensionales con las 
cuales obtendremos un sistema de ecuaciones mas compacto y con un signifi- 
cado fisico mas claro, de esta manera introduzcamos las nuevas variables: 

_ Q d _ 1 d 

3' d^'WoJt' ^ ' 

y las nuevas funciones adimensionales: 

-no -no -no t^c 

donde Hq es una constante que por conveniencia la escogeremos como 3Hq = 
k\ =^ \Ho\ = 1.66 X 10~^cm~\ notemos que en cosmologi'a Hq = 0.59 x 
10~^^cm~^ es la constante de Hubble cuyo inverso es una medida de la es- 
cala del universo [I/Hq = 1.69 x 10^^ cm), en nuestro caso tenemos que 
1/Ho ~ 6 km, magnitud razonable para nuestro sistema. 5*2 y S3 estan rela- 
cionados con la segunda y tercera componente del tensor cTq,/?. El nuevo 
tiempo adimensional (r) lo definimos a partir de la ecuacion fl4.1l) . notemos 
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que este tiempo puede tomar valores tanto positives como negatives, depen- 

diendo su signo del de Hq = i-i/— , para una mayor profundizacion en el 

V ^ 

significado del signo del tiempo adimensional ver el A. 

Comprobar que las variables antes definidas no tienen dimension es muy 
simple si recordamos que los coeficientes de la metrica (13 .Gp [A, B, C) son 
adimensionales, ademas de (I3.2ip y (I3.22p notamos que las variables O, S2 y 
S3 tienen unidades de inverso de longitud al igual que Hq. 

Si sustituimos las variables (14. 2 p en el sistema de ecuaciones (13.241) obte- 
nemos: 



^u,fj, L 



i2n - 53) (Fm - 2Tu,(3)p - m{Tp + r, 



(4.3a) 



52, r = -(3TM-3S2H, (4.3b) 

53, r = 2(3TM-3S^n, (4.3c) 

= fiTM-^^-\s2S,-h{'-\{Sl + Sl\ (4.3d) 
I3,r = 2(3{S-s-2n), (4.3e) 

la coma indica derivada respecto al tiempo adimensional (r), por ejemplo 
dS 

= — • Hemos cambiado la variable U por u, por (12.101) sabemos que 
dr 

U = U{/3, ^) ^ U^r = ^{^u,^l^^,T + ri/,/3/5,r) de donde podemos despejar /i^^. 
Ademas de (13.25P tiene que cumplirse que: 

- Si - Si - S2S3 + Sn^ = STu. (4.4) 

El sistema de ecuaciones diferenciales (14. 3 p lo resolveremos numericamente 
empleando el programa Maple. Sus soluciones nos daran el comportamiento 
en el tiempo del elemento de volumen del gas magnetizado de neutrones. De 
(14. 4p podemos despejar TC obteniendo dos raices, de las cuales escogeremos 
la negativa que garantiza (ver 14.51) la condicion de colapso. 



4.2 Soluciones numericas y discusion fisica 

4.2.1 Comportamiento de las funciones. 

Las soluciones que buscamos son las que nos dan el colapso del gas mag- 
netizado de neutrones que podria implicar el colapso de un volumen de la 
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ENs. Para garantizar la condicion de colapso basta con exigir en las ecua- 
ciones( H75i) que la expansion inicial sea negativa (7Y(0) < 0) pues su signo es 
el que nos indica si el volumen esta colapsando o expendiendose. De (14.21) y 
(I3.2ip se obtiene que el promedio del volumen local {V = ABC) lo podemos 
expresar como: 



Segiin el comportamiento de los coeficientes de la metrica podemos clasificar 
las singularidades como [T7] : 

1. tipo punto {A, B,C ^0), 

2. tipo cigarro (dos de los coeficientes metricos tienden a cero y el tercero 
se va a infinito), 

3. tipo barril (dos de los coeficientes metricos tienden a cero y el tercero 
tiende a un valor finito), 

4. tipo pancake (uno de los coeficientes metricos tienden a cero y los otros 
dos tienden a un valor finito). 

Estas denominaciones se refieren al cambio de forma del elemento de volumen. 
Fueron introducidas por Thorne (1967) en estudios sobre modelos Bianchi I, 
pero pueden aplicarse a casos generales. 

Para investigar el tipo de colapso a partir del comportamiento de los 
coeficientes espaciales de la metrica observemos que estos y las variables 
dinamicas antes definidas (7i + Si,i = 1, 2, 3) estan relacionados a traves de 
las siguientes expresiones: 



donde Aq, Bq, Cq son constantes, estas relaciones se obtienen de (13.211) y 



Para la solucion del sistema (14.31) utilizaremos diferentes condiciones ini- 
ciales tipicas para las ENs, por ejemplo: /i = 2 ^ p ~ lO^^g/cm^, Pq = 
10~^ — 10~^ para campos de entre lO^^G y lO^^G. Impondremos siempre 
la condicion de colapso del volumen 7i(0) < y tomaremos 5*2(0) = 0, ±1, 




(4.5) 




(4.6b) 
(4.6c) 




(E22]). 
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5*3(0) = 0, ±1 correspondientes a casos con deformacion inicial cero y defor- 
macion inicial en la direccion de los e]es y o z. 

La solucion de ?i(figura l4.ip para distintas condiciones iniciales nos mues- 
tra que 7i —00 independientemente de las condiciones iniciales. El campo 
magnetico tiende a aumentar manteniendose por debajo del campo critico. 
Su comportamiento para las condiciones iniciales estudiadas se puede apre- 
ciar en la figura 14.21 



■H 



- a=(s^(0)=-l, s^(0)=l, n=2, p(0)=0.01, w(0)=-2.46) 

- b=(s^(0)=-l, s^(0)=-l, ti=2, P(0)=0.001, W(0)=-2.59) 

- c=(s (0)=0, s (0)=0, 1^=2, P(0)=0.0001, H(0)=-2.39) 



— I — 

0.01 



— I — 

0.03 



• ■■■\ 

i !! 



Figura 4.1: Comportamiento de TL vs t para diferentes condiciones iniciales. 



De las soluciones del sistema de ecuaciones (14. 3 p obtenemos los graficos 
(14. 3p . (14. 4p cuyo comportamiento nos da la forma en que colapsa el elemento 
de volumen. 

Podemos observar como los Si + Ti se van a —00 por lo que los coeficientes 
espaciales de la metrica tienden a cero {A, B,C ^ 0) de donde podemos 
inferir que el volumen colapsa en una singularidad tipo punto. 
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p 



- a=(s^(0)=-l, s^(0)=l, 11=2, |3(0)=0.01, W(0)=-2.46) 

- b=(s^(0>-l, s^(0)=-l, n=2, p(0)=0.001, W(0)=-2.59) 
-c=(s (0)=0, s (0)=0, n=2, p(0)=0.0001, H(0)=-2.39) 



I » * 11 » « 



■ p ■ * * « ^« w * ^ ^ 

0.00 aa 0.02 0.03 0.04 O.OS 0.06 0.07 



' I I 

0.08 0.09 0.10 



Figura 4.2: El comportamiento del campo magnetico [f] = B/Bc) muestra una tendencia 
a incrementarse, manteniendose el valor del campo por debajo del campo crftico. 



■iiiiiii 



+ 



- a=(s^(0)=-l, Sj(0)=l, n=2, P{0)=0.01, ■H(0)=-2.46) 

- b=(s^(0)=-l, Sj(0)=-1, |a=2, |i{0)=0.001, H[0)--Z59) 
-c=(s (0)=0,s (0)=0,n=2, P(0)=0.0001, H(0)=-2.39) 



— I — I — I — 

0,02 0.03 



— I — I — I — 

0.04 0.05 



— I — I — I — I — I — 

0.07 0.08 0.09 0.10 



-ISO' 

+ 



- a=(s^(0)=-l, sJQ)=l, n=2, K0)=0,01, «(0)=-2.46) 

- b=(s^CO)=-l. 83(0)=-'. H=2, K0)=0.001, H(0)=-2.59) 

- c=(s (0)=0, s (0)=0, n=2, p(0)=0.0001, H(0)=-2.39) 



0.00 0.a 0.02 0,03 



—1 — I — I — I — I — I — I — I — I — 

0.04 0.05 0.06 0.07 008 



Figura 4.3: Comportamiento de {Si + Ti) y de (S'2 + TL) vs r. Podemos observar 
la tendencia a —00 ademas de distintos tiempos de colapso para diferentes condiciones 
iniciales 
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^ -150. 



SeiiHIIIKIIIIttllllllllli;..,.,., 

* * * • . 

• ■ t 

: !l 
: : I 



■ — a=(s^(0)=-l, s^(0)=l, n=2, P(0)=0.01, «(0)=-2.46) 
• — b=(s^(0)=-l, s^(0)=-l, ^=2, P(0)=0.001, H(0)=-2.59) 
c=(s (0)=0, s (0>0, |j,=2, P(0)=0.0001, H(0)=-2.39) 



-I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

0.00 001 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 



Figura 4.4: Grafico de (53 + Ji) vs r. Igualmcntc so observa que (^3 + 7Y) 



4.2.2 Espacio de fase. 



Si utilizamos la ecuacion (14.41) podemos reducir el sistema (14. 3 p a un sistema 
de ecuaciones en las variables 5*3, /5, H, quedando: 



1 



2(3{Ss-2n), 



(4.7a) 
(4.7b) 
(4.7c) 
(4.7d) 



notemos que la linica ecuacion que ha cambiado es la de Ti. 

Las trayectorias en la seccion {S^,[3,jj) del espacio de fase se muestran 
en la figura (14.51) . La evolucion del sistema queda determinada por el signo 
de Hq. Para r < =^ Hq = —^JnX/'i el sistema evoluciona hacia el atractor 
estable (punto marcado con la a), mientras que para t > Q ^ Hq = kX/3 
la evolucion es hacia una singularidad. Un estudio similar se realize para 
las demas secciones del espacio de fase obteniendose que las coordenadas del 
atractor son: (S'3 = 0, /3 = 0, /i = 1, 7i = 0). 
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5i 




Figura 4.5: Trayectorias en la seccion del espacio de fase (S'a, /?,//). Los puntos 
sin marca representan las condiciones iniciales, el punto "a" representa el atractor 



4.3 Conclusiones. 

Hemos estudiado la evolucion de un volumen local de un gas magnetizado de 
neutrones bajo condiciones que pueden ser las que existan en el interior de 
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una ENs. Para aproximarnos a la solucion de este problema hemos utilizado 
las ecuaciones de estado de diclio gas y las ecuaciones de Einstein-Maxwell 
junto a la ecuacion de conservacion de la energi'a, obteniendo de esta forma un 
sistema de ecuaciones diferenciales no lineales. El sistema obtenido lo hemos 
puesto en funcion de nuevas variables que permiten reducir el orden de las 
derivadas a la vez que dan una mejor interpretacion de los resultados. Para 
la solucion del sistema de ecuaciones diferenciales hemos introducido nuevas 
variables adimensionales, realizando los calculos numericos con el programa 
Maple. 

Los principales resultados obtenidos son: 

1. La dinamica del volumen magnetizado de neutrones presenta una solucion 
singular que equivale a un colapso de la materia. Este resultado se 
obtiene para todas las condiciones iniciales estudiadas. Con ello res- 
pondemos la pregunta de que si los efectos de la RG al ser tomados en 
cuenta favorecerian o no el colapso magnetico descrito en [2]. Como 
podemos ver la gravitacion hace que el colapso se manifieste y ocurra 
en forma de punto. 

2. A diferencia del estudio dinamico realizado para el gas de electrones 
magnetizado [3] el tipo de singularidad que aparece en este caso es de 
tipo punto y no cigarro. Este resultado sugiere que los efectos de la 
gravedad para un gas de neutrones magnetizado son mas marcados que 
para el gas de electrones, cambiando la forma de la singularidad. 

3. El espacio de fase muestra que el sistema evoluciona, para Hq < 0, a 
un punto de equilibrio, es decir una configuracion estable. 

4. El campo magnetico tiende a aumentar, manteniendo valores por de- 
bajo de campo critico. 

4.3.1 Direcciones del trabajo future 

Resultan'a interesante estudiar un volume de gas de neutrones, electrones 
y protones que obedezcan el equilibrio /? y la neutralidad de carga, lo que 
describirfa de manera mas realista las condiciones en el interior de una ENs. 

Otros espacio-tiempo como Bianchi V, VII 6 IX permitirian tomar en 
cuenta mas grados de libertad, acercandose mas a la descripcion de lo que le 
ocurre a un volumen local de una ENs. 
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Anexo A 



Significado del tiempo 
adimensional (r). 



El tiempo adimensional (r) se define a partir de la ecuacion fl4.1l) : 

> r = Hot, (A.l) 



d _ 1 d 
dr Hq dt 



en la cual podemos notar que el signo de r depende del signo de t y del 
signo de Hq. Para comprender el sentido fi'sico del tiempo adimensional r 
observemos que de (14.21) y (13.211) obtenemos: 



HHq = H 



1 A B C 1 d 
S^A^B^C^ ~ 3di 



HV), 



V =ABC^ — 




(A.2a) 
(A.2b) 



y despejando H de (14. 4p tenemos: 



n = ±-^/3Tu + si + si 



SoS: 



2^3) 



(A.3) 



Si fijamos t ^ 0, de ( ]A.2|) tenemos la siguiente interpretacion fisica: 



1 = 3 J Hdt 





I>O^V>Vq^ expansion, 
/<O^V<V'o^ colapso 
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Las posibles combinaciones de signo se dan en el cuadro A, en el cual 
aparecen resaltados los convenios que hemos escogido en esta tesis. 



Casos 


Expansion 


Volumen elemental 


T 


H>0 


n>Oj Ho>0 


Expansion 


r > 




HkOj Ho<0 


V>Vo 


r < 


H <0 


n<OjHo>0 


Colapso 


r > 




n>OY Ho<0 


V<Vo 


r < 


H = 


n = o 


V = Vo 





Tabla A: Posibles combinaciones de signos entre Ti, Hq j t para t > 0. 
Notemos que el sistema comienza a evolucionar a partir de t = r = 0, para 
r > el sistema colapsa y para r < el sistema tiende hacia un estado 
diluido como se puede ver en el grafico (14. 5p . 
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